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We argue that the neutrino advance of time observed in MINOS and OPERA experiments can 
be explained in the framework of the standard relativistic quantum theory as a manifestation of the 
large effective transverse size of the eigenmass neutrino wavepackets. 
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INTRODUCTION 

Recently, OPERA Collaboration reported that neutri- 
nos from CERN arrive to the Gran Sasso Underground 
Laboratory by (60.7 ± 6.9stat. ± 7.4sys.) ns earlier than 
expected for almost massless particles [l|. MINOS Col- 
laboration also observed f^] an earlier arrival of neutrinos 
from FNAL to the Soudan Underground Laboratory by 
(126 ± 32stat. ± 64sys.) ns (68% C.L.). The experiments 
have similar distances of about 730 km between the neu- 
trino production and detection regions but different mean 
neutrino energies (17 and 3 GeV for, respectively, CNGS 
and NuMI beams) and different neutrino flavor compo- 
sitions of the beams. These remarkable results, being 
interpreted in terms of neutrino velocity Vi,, suggest a 
superluminal motion of neutrinos with 



(5.1 ±2.9) X 10"^ (MINOS), 
(2.48 ± 0.41) X 10"^ (OPERA). 



At the first blush this interpretation breaks the relativity 
principle - one of the basis of modern physics. We will 
however argue that an earlier arrival of neutrinos could 
be understood without any violation of relativity, causal- 
ity and other fundamental physical concepts and is just 
a manifestation of quantum nature of neutrino. Namely, 
we will try to demonstrate that the observed effect can 
be at least qualitatively explained by taking into account 
that the quantum states of neutrinos with definite (small) 
masses are described by the relativistic wavepackets hav- 
ing a finite and in fact very large effective transversal size. 
Necessity of a wavepacket description of neutrino prop- 
agation in vacuum and matter is now well understood 
(while not yet commonly accepted) in the theory of neu- 
trino fiavor transitions ( "oscillations" ) based on quantum 
mechanics or quantum field theory. 



NEUTRINO WAVEPACKET 

Any wavepacket can be conventionally characterized 
by a most probable 4-momentum, 4-coordinate, and a 
set of parameters governing the shape of the packet in 



the phase space. Apparently, a spherically symmetric 
wavepacket with an effective spatial "size" ax and mo- 
mentum "width" cTp ^ ^/(Jx in its rest-frame becomes 
asymmetrical if it is boosted with a Lorentz factor F ^ 1. 
The wavepacket spatial size in the boost direction shrinks 
as Gx/r remaining unchanged in the transverse plane. 
The momentum width increases in the boost direction 
as CTp-T remaining the same in the transverse plane. In 
our previous paper [3] we developed a covariant field- 
theoretical approach to neutrino oscillations which oper- 
ates with the relativistic wavepackets describing initial 
and final states of particles involved into the production 
and detection of neutrino. The neutrino in this approach 
is described as a virtual mass eigenfield travelling be- 
tween the macroscopically separated production and de- 
tection vertices of Feynman graphs. Thus we make no 
any assumption about its wavefunction. Instead, within 
our formalism we compute the neutrino wavefunction 
which turns out to be a wavepacket with spatial and 
momentum widths defined and functionally dependent 
on those of the particles involved into the neutrino pro- 
duction and detection subprocesses. Explicitly, up to a 
coordinate independent spinor factor, the effective (out- 
going) neutrino wavefunction reads 



c)-<T^r„^(x||-V„Xo)'-<T,^ 



(1) 



Here xq is the time, X|| and x_l are, respectively, the lon- 
gitudinal and transverse (relative to the mean velocity 
vector -Vy = p^/E^) spatial coordinates of the geometric 
center of the neutrino packet {xq and X|| -I- xj_ form a 
4- vector); Ei, = \/pl + ml, and nii, is the mass of the 
neutrino mass eigenfield. In the most general case, the 
"spread" ai, is a Lorentz invariant function of the most 
probable 3- momenta p^, masses m^t, and momentum 
spreads — const of all external in and out particles 
x: involved into the neutrino production-detection pro- 
cess which are described as asymptotically free relativis- 
tic wavepackets. It is shown in Ref. that the center 
of any external wavepacket moves in the mean along the 
classical trajectory (x;^^) = Sc^^ + v^x^ conserving energy, 
momentum and effective volume (oc 1/crJr); under certain 
conditions the packets remain stable (nondiffluent) dur- 
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ing the times much longer than their mean lifetimes (in 
case of unstable particles) or the mean time between the 
two successive collisions in the relevant ensemble (in case 
of stable particles) . 

Considering that the two-body decays of pions and 
kaons are the main processes of neutrino production in 
the MINOS and OPERA experiments, we can neglect the 
contributions into coming from the particles, interact- 
ing with neutrinos in the detector (reasonably assuming 
that their 4- momentum spreads are much larger than CTtt , 
gk and cr^). With this simplification, it can be proved 
that 



AB = AD = i 



Neutrino 
\ wavepacket 





TT or K. 



Then from the above-mentioned conditions of stability 
for the meson and muon wavepackets it follows that 
must satisfy the following conditions 



al < ml 



r„ 



where Ta = 1 /Ta and — 1 /r^ are the full decay widths 
of the meson a and muon. Considering that for any know 
meson m^/Tf^ ^ ma/Ta, we conclude that 



al/ml < r^/m^ w 2.8 x 10 



-18 



(2) 



Therefore the neutrino momentum uncertainty is fantas- 
tically small [g^. From ([2]) one can immediately derive 
the lower bounds for the effective spatial dimensions of 
the neutrino wavepacket: 



= ^ » 10-^ 



0.1 eV 

nil, 
10 GeV 



km, 
0.1 eV 



^m, 



So the neutrino wavepacket appears as a huge but su- 
perfine disk of microscopic (energy dependent) thickness 
in longitudinal direction, comparable with the thickness 
of a soap-bubble skin, and macroscopically large (energy 
independent) diameter in the transverse plane . 

This is a key point in interpretation of the experiments 
with earlier arrival of neutrino signal. 



QUALITATIVE ESTIMATIONS 

In fact, neutrinos produced at accelerators arrive to the 
detector site widely distributed across the plane trans- 
verse to the beam axis. Those neutrinos which were mis- 
aligned with the neutrino detector nevertheless do have a 
chance to interact within the detector due to macroscopi- 
cally large transverse size of its wavefunction. Moreover, 
its interaction probability very weakly depends on the 



FIG. 1. Neutrinos are emitted from the "Somce" and are 
registered in the "Detector". The centers of the neutrino 
wavepackets will arrive at the points B and D simultane- 
ously, while the signal from the neutrino wavepacket (shown 
as an extremely oblate spheroid) which moves under the angle 
6 = ZBAC to the beam axis will arrive earlier since DE > 0. 
Neutrino velocity vector v^ lies in the plane of the figure. 
Proportions do not conform to reality. 



misalignment distance r = BC (see Fig. [T|) if it is small 
compared to d±. As is seen from Fig. [TJ the misaligned 
neutrinos will interact systematically earlier than those 
moving along the beam axis, due the huge transverse 
width of their wavefunctions. The school-level planime- 
try suggests that the advancing time is given by 



(5i = L(l/cos6i-l) «rV(2L). 



(3) 



Here we assume that (i) 1 — Vi, ^ 1, (ii) the neutrino 
wavepacket effective width is much larger than the de- 
tector dimensions, and (iii) 1. Substituting num- 
bers into ([3]) one obtains Fig. [2] from which quantitative 
estimates for the time advance as a function of r could 
be drawn. For instance, a neutrino packet which moves 
3 km away from the detector will come earlier by about 
20 ns than that moving directly to the detector. 

What is the probability to find a neutrino at a distance 
r from the beam axis? This could be estimated taking 
into account that neutrino production is dominated by 
two-particle decays of pions and kaons. The angular dis- 
tribution of massless neutrinos from these decays is well 
known: 



dl 
dh 



1 



1 



47r(l 



7r(l + r, 



(4) 



Here 9 is the angle between the momenta of the meson a 
and neutrino (0 < < tt), is the meson velocity, and 
Ta = (1 — w^)^^^^ = Ea/rHa- The second approximate 
equality in Eq. (j4]) holds for small angles and relativistic 
meson energies {9 ^ 1, AF^ 3> 1). In the latter case, 
the main contribution to the neutrino event rate comes 
from the narrow cone 9 < l/Ta- Considering that the 
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r (m) 

FIG. 2. Advance 5t as a function of r. 



mean neutrino energy, i?^, from the muonic decay of a 
meson with energy Ea is = raEi''\ where eI"^ = 
(m^ — m^)/(2ma) is the neutrino energy in the rest frame 
of the particle a, the characteristic angle can be defined 
as e^a) = Ei'^^/E,. 

In the case of OPERA, one can (very) roughly estimate 
the characteristic angles for the "low-energy" (LE) range 
{E^ < 20 GeV, E^ ^13.9 GeV) and "high-energy" (HE) 
range {Ei, > 20 GeV, E^ 42.9 GeV), assuming that the 
main neutrino sources in these ranges are, respectively, 
7r^2 and i4r^2 decays: 

^LE > V) = 2.1 X 10-3, 
She < 0iK) = 5.5 x 10-^. 

This provides us with an ordcr-of-magnitude estimate of 
the mean values of r and advancing times St: 

rhE > 1-7 km, the < H km; 
SthE ^ 5.6 ns, StuE ^ 36.7 ns. 

Since the LE and HE ranges contribute almost equally 
to the CNGS t'^ beam, there must be a definite trend 
towards earlier neutrino arrival to OPERA with approx- 
imately 21 ns mean time-shift and a "tail" or, better to 
say, "fore" of the same order coming from the "edges" of 
the CNGS beam. 

Similar estimation for the low-energy NuMI beam at 
Fermilab producing neutrinos for the MINOS experiment 
can be done with a better accuracy, since the 7r^2 decay 
is here the dominant source of neutrinos and the radial 
distribution of the beam is expected to be very flat. So, 
by using E^ = 3 GeV we obtain 

r « 36.2 km, 6t w 120.7 ns. (5) 

The latter number is in surprisingly good agreement with 
the MINOS observation. Obviously, MINOS should ob- 
serve at the average a much earlier arrival of neutrinos, 
in comparison with OPERA, because of the lower mean 



neutrino energy which corresponds to a wider transverse 
beam distribution and hence to a larger input from the 
misaligned neutrinos. 

NUMERICAL RESULTS 

Let us reevaluate the estimates given above with a 
somewhat detailed but still simplified calculation. In par- 
ticular, we could profit from the simulation of expected 
radial distribution of i/^ charged current (CC) events per- 
formed by the OPERA Collaboration This distribu- 
tion {Pccif)) which we digitalized for our purposes is 
displayed in Fig.[3l Being dominated by the 7r^2 and Kf^2 
decays, the transverse beam size at Gran Sasso is of the 
order of kilometres and the full width at half maximum 
of the distribution is about 2.8 km. 
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FWHM « 2.8 km 
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FIG. 3. Probability of neutrino charged current interactions 
expected in OPERA as function of r. 

The distribution Pcc('') transformed (with help of 
Eq. ([U) into the 6t distribution as 

p rPcc {r jSt)) 

Pcc[Ot) = -Joe 

/ drrPcc (r) 

JO 

is shown in Fig. ID Its average {6t) is about 20 ns with 
similar variance and with the tail extending up to about 
100 ns. Figure [5] shows the integral distribution 

Pcc{<St) = / dt'Pccit'). 

Examination of Fig. [5] suggests that all CC events roughly 
equally populate the following intervals in St: (0, 20) ns, 
(20,45) ns, and (45, 100) ns. 
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FIG. 5. Pcc(< St) distribution expected in OPERA. 



FIG. 6. Top panel: time probability density function for the 
first beam extraction taken as an example, expected by the 
OPERA Collaboration 1] and our calculation. In both cases, 
the systematic "instrumental" shift is not shown since it does 
not change the shape of the curves. Bottom left and right 
panels: zooms of the top panel for the left and right fronts of 
the signal, respectively. 



In the likelihood fit of the time distribution performed 
by the OPERA Collaboration the fronts of the time dis- 
tribution statistically play the major role. Therefore, the 
time distortion evaluated in the present work seem to 
explain the OPERA observations without any model pa- 
rameter and without introducing superluminal neutrinos 
or other exotics. 



Finally, we compute the expected time distribution in 
OPERA, g{t), as a convolution of the probability den- 
sity function of arrival time f(t) shown in top panel of 
Fig. [6] as solid line, taking into account an earlier arrival 
of neutrino signal as follows: 

/•OO 

/ f {t + St(r)) Pcc{r)rdr 
m = ^-^ p . 

/ Pcc{r)rdr 
Jo 

The resulting curve g{t) is displayed superimposed in the 
top panel of Fig. |6] by dashed line. As is seen, on the 
average, it is shifted to the left by about 20 ns. However, 
and this is even more important, the left front of the 
signal is shifted by a larger amount as it accumulates 
the advance effect from the total f{t) distribution. One 
could verify that the solid and dashed lines shown in the 
bottom panel of Fig. |6] for the left front of the OPERA 
time distribution is shifted by about 50 ns (up to 60 ns on 
the tail) to the left, while the right front is shifted only by 
20-25 ns. In other words, the impact of the misaligned 
neutrinos is predicted to be asymmetric in time. 



CONCLUSIONS 

Large transverse size of the neutrino wavepacket and 
uncollimated beam of neutrinos seem to explain the ear- 
lier arrival of the neutrino signal in OPERA and MINOS. 
The neutrino signal is estimated to arrive in advance by 
about 20 ns in the mean (with a similar variance) for 
OPERA and by about 120 ns for MINOS. In the case 
of the OPERA experiment only this effect essentially re- 
duces the statistical significance of its observation. More- 
over, we have evaluated the expected time distribution of 
the neutrino arrival in OPERA and obtained that the left 
and right fronts are shifted to the left by about 50-60 ns 
and 20-25 ns, respectively. This probably explains the 
observed anomaly all-in-all without any exotic hypoth- 
esis, like Lorentz violation and so on. Let us underline 
that in our calculations we do not use any adjustable 
parameter. In the case of the MINOS experiment there 
is also a surprisingly good agreement between our ex- 
pectation ^ and experimental result. Therefore, we ar- 
gue that observations of superluminal neutrinos by the 
OPERA and MINOS experiments can be treated as a 
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manifestation of the huge transverse size of the neutrino [1] 
wavefunction. This kind of effects could be investigated 
in the future experiments (in particular, in the off-axis I^l 
neutrino experiments) with more details in order to prove 
or disprove our explanation. 

Let us note that one should not expect an increase in 
the number of neutrino induced events due to the mis- 
aligned neutrino interactions because this effect will be [4] 
compensated by the corresponding decrease of the num- 
ber of aligned neutrinos. 

Let us briefly discuss the situation with the observed 
(anti)neutrino signal from SN1987A. A proper treatment jg-j 
of these neutrinos should take care about the dispersion 
of the neutrino wavepackets at astronomical distances. 
Deliberately neglecting the dispersion, it appears that 
any terrestrial detector is sensitive only to the aligned 
neutrinos, since the misaligned neutrinos will have nee- 
ligible impact due to the smallness of their wavepacket ^ ' 
transverse size relative to the astrophysical scale of about 
50 kps. Therefore, no advance signal should be expected. 
However this problem is not so simple and needs in a 
more detailed theoretical analysis. 
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In fact expression ([T]) has a limited range of applicabil- 
ity which is however deliberately comprises all terrestrial 
neutrino experiments. 

By the way, Eq. explains a success of the naive 
quantum-mechanical approach to neutrino oscillations 
which is explicitly based on the assumption that the mas- 
sive neutrinos have definite 3-momenta and thus can be 
described by plain waves (infinitely large mathematical 
artifacts) . 

Note that the acceptable effective transverse dimensions 
of "normal" particle states (or, that is the same, the ef- 
fective diameters of the packets in the intrinsic frame of 
reference) are not so huge as for neutrinos. For example, 
the limiting dimensions (minimal effective diameters of the 
wavepackets, necessary to keep them stable during the life 
of the particles) should be of the order of 10~^, 10"^, and 
10~* cm for TT*, fi^, and r*, respectively. These dimen- 
sions are probably hard to measure in the current particle 
scattering experiments, since the wavepackets are reduced 
due to any measurement (interaction). Neutrinos provide 
the only (as yet) tool for such measurements, which are 
believed to be already realized by the MINOS and OPERA 
experiments. 
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